Energy Relaxation in one-dimensional polariton condensates 
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We study the kinetics of polariton condensation accounting for the condensation process as well as 
the energy relaxation of condensed polaritons due to their scattering with phonons and excitons. By 
assuming a Boltzmann kinetic description of the scattering process, we show that intra-condensate 
relaxation can be accounted for by an additional time-dependent term in the Gross-Pitaevskii equa- 
tion. As an example, we apply the formalism to the experimental results recently obtained in 
polariton microwires [E. Wertz, et al., Nature Phys. 6, 860 (2010)±]. In the presence of a local 
non-resonant optical pump, a dynamic balance between spatially dependent relaxation and particle 
loss develops and excites a series of modes, roughly equally spaced in energy. Upon comparison, 
excellent agreement is found with the experimental data. 
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I. INTRODUCTION 

Researchers now have two decades of experience in 
the study of semiconductor microcavities^ and exciton- 
polaritons, the mixed quasiparticles that arise from 
strong light-matter coupling. Inspired by the early pre- 
diction of Bosc-Einstcin condensation of excitons^r— , im- 
pressive progress in cold-atom condensation^, and mo- 
tivated by the prediction of ultra-low threshold polari- 
ton lasing2£, many experiments have been dedicated to 
the realization of a Bose-Einstein condensate of polari- 
tons£~— . As mixed quasiparticles, polaritons are ex- 
pected to interact with each other (due to underlying 
exciton-exciton interactions) whilst maintaining a light 
effective mass (due to the trait of cavity photons). As 
bosons, polaritons should relax into the ground state of 
the system, producing a transition to a phase coherent 
state capable of a laser-like emission without the need 
for population inversion^— . This transition is the non- 
equilibrium analog of Bose-Einstein condensation and 
can occur at high temperature due to the light effective 
polariton mass. Although research studies usually use 
an optical excitation of polaritons, progress in electrical 
injection sets a direct path for practical applications in 
the form of polariton-bascd light-emitting diodes^—. A 
growing trend towards studying highly confined struc- 
tures, in which polaritons are confined into zero or one 
dimensional structure o 14 ' 20 " — , opens the way to the re- 
alisation of microstructured polariton devices. 

A theoretical description of the condensation of mi- 
crocavity polaritons is a challenging 2 ^— task, especially 
when one wishes to account for spatial dependence. Since 
the polariton lifetime is shorter than the relaxation time 
(in today's samples), the system always lies in a state 
out of thermal equilibrium, with no well defined chemi- 
cal potential or temperature. For this reason one requires 
a kinetic model that accounts for the excitation, scatter- 
ing, relaxation and decay processes occurring in the sys- 
tem. Semiclassical Boltzmann rate equations^ can be 
used to describe the thermal relaxation processes that 
lead to condensation, however they completely lack of a 



description of the collective excitations that characterize 
a quantum fluid. To correctly describe these properties 
at thermal equilibrium, symmetry-breaking theories are 
instead necessary. These approaches model the conden- 
sate as a classical field, in terms of the Gross-Pitaevskii 
equation, possibly coupled to a quantum model for the 
excited states, in analogy to the theory of atomic conden- 
sate a 24 ' 26 ' 29 . Intuition suggests that we can use a mix of 
the two theories in which higher energy excitons are de- 
scribed by Boltzmann-like equations and condensed po- 
laritons by Gross-Pitaevskii-likc equations. Such a model 
was constructed by Wouters and Carusotto 3 ^ to describe 
polariton condensates, whilst similar approaches have 
been developed in atomic condensates 3 ^ - — . 

Still, in experiments condensation can take place into 
a state above the ground stat o 34 ' 35 and in planar mi- 
crocavities condensed polariton relaxation was evidenced 
down a potential gradient 1 ^. In periodic structures, con- 
densation was observed into both in-phase and antiphase 
Bloch states and the further relaxation into the lower en- 
ergy in-phase state demonstrated with increasing pump 
intensitj*ii. Particularly stunning relaxation effects have 
appeared in recent experiments on polariton microwires 1 -. 
Here condensation was observed into a series of energy 
modes that appear roughly equally spaced in energy. Sur- 
prisingly, the energy spacing is too large to be attributed 
to quantization of the center-of-mass motion (the mi- 
crowire is too long). To describe these effects theoret- 
ically it is essential to have a model that accounts for 
both the energy relaxation of excitons into condensed 
polaritons and relaxation of polaritons within the con- 
densed fraction of the polariton gas (polaritons can relax 
their energy by scattering with phonons or hot excitons). 
While the first order process is included in most kinetic 
theories of polariton condensation, the second has never 
been theoretically addressed, to our knowledge. 

In this work we derive a generalized version of the 
Gross-Pitaevskii equation (GPE), that describes intra- 
condensate relaxation induced by the coupling with a 
reservoir. The starting point of the derivation is the 
Boltzmann collision term for the condensate density. As 
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a result, the GPE acquires an additional time-dependent 
rate that models the relaxation process. As an example, 
we apply the formalism to the recent experiment carried 
out by Wertz et al. on polariton microwiresi The in- 
coherent optical pump is modeled using a classical rate 
equation coupled to the generalized GPE. The optical 
pump itself introduces a spatially structured potential 
(an effect also shown in experiments using resonant ex- 
citation^). Only a full account of the energy relaxation 
within the condensate and of the spatial dynamics of po- 
laritons interacting with this potential can reproduce the 
various modal patterns measured experimentally. 

The paper is organized as follows. In Sec. [HI we derive 
the term in the GPE that describes the relaxation. Sim- 
ulations of relaxation in nonresonantly excited ID wires 
are presented in Sec. IIIII Conclusions are presented in 
Sec. EH 



II. DESCRIPTION OF THE RELAXATION 

The degeneracy of the Bose gas implies that a good 
account of the condensate kinetics is already given by a 
mean field description, formulated in terms of a classical 
field ip(x,t). This approach is common to most recent 
theoretical works^— that describe most existing exper- 
imental data obtained above the condensation threshold. 
More advanced quantum kinetic theories are only needed 
to model the details of phase fluctuations 4 ^—. 

A kinetic model formulated in terms of transitions be- 
tween states characterized by their frequencies, involves 
a coarse graining of the time. This is analogous to the 
coarse graining of space used in the semiclassical Boltz- 
mann equation 4 ^ to describe the density as a function 
of position and momentum while fulfilling the constraint 
imposed by the uncertainty principle. We will use the 
following definition of the time dependent spectrum of 
the bose field ip(x,t): 



1 

T 



ip(x,t')dt', 



(1) 



t-T 



where T is the coarse graining time step. When using 
this formalism to model a steady-state spectrum, it is 
important that h/T is chosen smaller than the energy 
spacing of the features to be resolved. Then, the physical 
result should not be affected by the precise value of T. 
The inverse transform reads 



(2) 



Spontaneous scattering is a quantum feature that can- 
not be described within a classical field model such 
as the GPE and has only been treated in a handful 
of theorie s 46 ' 47 ' 49 . Neglecting spontaneous scattering 
might limit the accuracy in modeling the kinetics of low- 
occupation states. It is however a particularly well suited 



assumption when modeling intra-condensatc kinetics. In- 
deed, due to the high occupation of both initial and fi- 
nal states in a scattering process, stimulated scattering is 
largely dominant. We therefore restrict to the stimulated 
relaxation processes. At high densities, relevant for po- 
lariton condensates, these processes are mainly polariton- 
exciton scattering^ rather than polariton-phonon scat- 
tering^. The in-scattering term at frequency w reads 



dn LU {x,t) = n u {x,t) ^^r{uj,uj')n u >{x,t) dt, (3) 



where r{u!,u>') is the net scattering rate from the mode 
at frequency u>' to the mode at u. Due to particle num- 
ber conservation, the relaxation rate obeys the relation 
r{uj,uj') = — r(uj',Lu). The first term in the expansion 
of the rate as a function of the frequency difference is 
therefore 



(4) 



where the relaxation constant k has the dimension of an 
inverse density. Here, we don't use a microscopic model 
of the scattering process. Rather, we keep free 
parameter that we adjust to fit the experimental data. 
The resulting value is consistent with the scattering on 
hot excitons as discussed by Porras, et ali££. 

Under the change of the density ([3]), the wave function 
varies as 



4> u {x,t) + dip u (x,t) = 



'n u {x,t) + dn^{x,t) 



n u {x,t) 



(5) 

where we used the fact that a stimulated relaxation pro- 
cess does not change the phase of ip. Expanding Eq. ([5]) 
to first order in dn^, substituting Eqs. ([3]) and ((4]) and 
using the inverse transform ([2]), one obtains for the re- 
laxation dynamics of the wave function 



dip{x,t) nn(x,t) 
dt = 2 



(i(x,t) 



id 
di 



i>(x,t). (6) 



In the derivation, boundary terms in the partial integra- 
tion that scale as 1/T were neglected, because we con- 
sider the limit of a large coarse graining time where 1/T 
is smaller than any other relevant frequency scale. In 
Eq. ([6]), n and /2 are defined as 
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\^{x,t')\ 2 dt' 
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id 



dt 



(7) 
.(8) 



The first expression defines the density averaged over 
the coarse graining interval T. Analogously, hp, can 
be suggestively read as an averaged chemical potential. 
We stress however that no actual chemical potential can 
be defined for this non-equilibrium system, and Eq. J5J) 
draws its only justification from the formal derivation 
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presented above. The real part is taken in Eq. (JSJ, be- 
cause it is readily shown with the Madelung transforma- 
tion if) = y/n e lS , that the imaginary part of the integral 
in Eq. ([5]) scales as 1/T. Equation ^ gives the term 
that was sought to describe the stimulated relaxation of 
a classical Bosc field due to local scattering with phonon 
like particles and is the central result of this work. 

The right hand side in Eq. © resembles the fre- 
quency dependent amplification term that was intro- 
duced in Ref. [Hi] to describe an energy dependent gain 
mechanism due to scattering from a reservoir: dip/dt = 
(P/fijf)[Qir — id/dt]ip. Two main differences should be 
emphasized. First, here the relaxation term is propor- 
tional to the polariton density n, whereas in Ref. [5l|, the 
relaxation is proportional to the gain from the reservoir 
P. The second important difference is that the gain cut- 
off frequency VLk is replaced by the average frequency 
p,. As a consequence, particle loss and gain balance each 
other and there is no net gain in Eq. ([5]). The modes 
with frequency u> > fl experience loss, where the ones 
with uj < ft are amplified. Energy absorption from the 
reservoir into the condensate is thus not included in our 
model. The relaxation © is due to the interaction with 
an environment that is effectively at zero temperature. 



III. RELAXATION IN ID WIRES 

As an example of application, we consider a ID mi- 
crowiroi under non-resonant optical excitation. In most 
situations in which an incoherent optical pump is present, 
a good description of the condensation kinetics is given by 
coupling the GPE to a kinetic equation for the reservoir 
population; 3 — To this purpose, we introduce a complex 
condensate field, ip(x,t), coupled to a reservoir popula- 
tion, N(x,t). Neglecting the polarization degree of free- 
dom, the evolution of condensed polaritons is given by 
the generalized GPE for a Bose gas with losses and gain 
from a reservoir, supplemented with the dissipation term 
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+V(x,t)+g\il>{x,t)\' 



hnn(x, t) 



i^{x,t) + — 



iP(x,t), (9) 



where the kinetic energy eigenvalues of E^p are: 

E LP (k) = Ecik)+ 2 Exik) -y(Ec(k)-E x (k)f + ^. 

(10) 

Here Ec{k) and Ex(k) represent the cavity photon and 
exciton dispersion relations, respectively, and is the 
exciton-photon coupling constant. Polaritons scatter 
from the reservoir into the condensate at a rate Rr and 
radiatively decay at a rate Tc- V(x,t) = hRp_N(x,t) + 



HGP(x) is an effective potential experienced by polari- 
tons, caused by interactions between polaritons with the 
reservoir and also the pump field, P{x), described with 
a strength G. The large potential experienced by polari- 
tons at the ends of the wire is accounted for through the 
introduction of boundary conditions ip(x = ±L) = 0. 

The dynamics of the reservoir field, representing opti- 
cally injected hot excitons injected, is described by: 



dN(x,t) 
dt 



= - (r R + PR R \ip{x, i)| 2 ) N(x, t) + P(x) 



(11) 

where Tr is the decay rate of reservoir polaritons and 
the parameter (3 < 1 quantifies the backaction of the 
condensate on the reservoir—. Such a parameter is nec- 
essary due to the simplified model of a single reservoir 
used for the exciton states; if /3 = 1, then unphysical 
instabilities may develop in the solution of the coupled 
Eqs. © and (fTl) 46 . The reason is that for /3 = 1 the gain 
saturation mechanism is overestimated, given the simpli- 
fied description of the reservoir in terms of Eq. (jlip . At 
the same time, the assumption ft = would make the 
condensate density diverge, because the gain saturation 
mechanism is needed to determine the steady state con- 
densate density at a given pump intensity. Ultimately, 
a microscopic description of the reservoir made with a 
continuous spectrum of modes, would in principle make 
the introduction of this phenomcnological parameter un- 
necessary. The choice /3 < 1 is therefore obliged, given 
the reservoir model (fTTj) . Finally, note that the value 
of f3 makes no qualitative difference to the spectral and 
spatial condensate pattern (implying only an overall de- 
pendence of the condensate density) so long as it is small 
enough to prevent instabilities. 

We have solved Eqs. ([9lfTT]) with parameters^ corre- 
sponding to the experiment of Werta^ et al., who consid- 
ered two cases: excitation at the wire center and excita- 
tion near one edge. In both cases the time evolution is 
calculated over a long (4000ps) period during which the 
system reaches a steady state. Energy distributions are 
then calculated by Fourier transforming the fields. 



A. Excitation in wire center 

The energy distributions with excitation at the wire 
center are shown in reciprocal (left) and real space 
(right) in Fig. Q] The presence of the pump introduces 
excitons into the reservoir and also creates the spa- 
tially varying potential, V(x). Together with a slight 
blueshift caused by polariton-polariton interactions, this 
potential determines the energy at which a polariton 
condensate forms due to the scattering of particles from 
the reservoir. The condensed particles arc created with 
zero in-plane wavevector, in the center of the pumping 
spot, and are pushed away by the potential V(x) so 
that they fill the wire. As the pump intensity increases, 
a clear blueshift of the condensed polaritons is visible 
and stimulation of the energy relaxation takes over. 



Propagating polaritons form a spectral pattern with 
several modes, roughly equally spaced, in agreement 
with the experimental results^. 
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FIG. 1. (color online) For excitation in the wire center, the 
energy distribution of condensed polaritons in real (right) and 
momentum (left) space is shown for increasing pump powers. 
The white curves show the bare polariton dispersion given by 
Eq. (|10[) in momentum space (left) and the effective polari- 
ton potential, V(x), in real space (right). In each plot the 
intensity is scaled to the maximum intensity and plotted with 
a logarithmic color scale. The data were convoluted with a 
small spectral linewidth (8 — 0.02meV) to account for the 
finite spectral resolution in the experiment. 

The energy relaxation introduces a gain mechanism for 
the lower energy polariton states and one can imagine the 
relaxation of polaritons as they travel down the pump- 
induced potential, as illustrated in Fig. [5^. The rate of re- 
laxation between any two modes is an increasing function 
of their energy difference, AE, to lowest order according 
to Eq. as illustrated by the blue (dark-grey) line in 
Fig. [2J). For the formation of the condensate to occur in 
a lower energy mode, in the steady state, the gain must 
exactly balance the losses in the process. In other words, 
the incoming rate due to this relaxation must balance the 
outgoing rate, which is mostly due to radiative decay, il- 
lustrated by the green (light-grey) line in Fig. (2b. The 
balance occurs at the energy AEq and relaxation is only 
expected to appear between modes separated at this en- 
ergy. The further outward propagation of polaritons can 
allow a cascade effect in which polaritons relax in units 
of AEq into a series of levels (Fig. [2Jt) . As the rate is pro- 
portional to Kfi(x, t) and to the energy step, the modes 
are roughly equally spaced in energy if we assume that 
the density in the initial mode is approximately the same 
for all relaxation steps. The data by Wertz et al. support 
this. 
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FIG. 2. (color online) a) Illustration of quantized energy re- 
laxation as polaritons propagate down the pump-induced po- 
tential, V(x). b) Sketch of the incoming (proportional to 
Kfi{x,i)) and outgoing rates (r) for the lower energy mode in 
a relaxation step, as a function of the energy difference with 
a higher energy mode. 

To confirm the necessary role of energy relaxation in 
the generation of a series of condensate modes at dis- 
crete energies, we repeated our calculations setting k = 0. 
Fig. El shows that in this case all polaritons condense into 
a single energy mode determined mostly by V(x = 0). 
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FIG. 3. (color online) Same as the bottom panels of Fig. Q] 
with k = 0. 



B. Excitation near wire edge 

The case of excitation near the wire edge is shown in 
Fig. [4] The same phenomenology as in the case of exci- 
tation at the wire center can be identified; condensation 
at a blueshifted energy and energy relaxation at higher 
intensities. However, one can now observe the striking 
appearance of harmonic modes due to a potential trap 
formed between the wire edge and the pump-induced po- 
tential, V(x). This feature is in close agreement with 
the experimental results^. The positioning of the optical 
pump near the wire edge introduces a strong asymme- 
try between relaxation on the two sides of the pump. It 
is observed theoretically and experimentally that relax- 
ation is more effective when polaritons interact with the 
wire edge. 

We note that our model is limited in its quantitative 
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FIG. 4. (color online) Same as in Fig. [T] for excitation near 
the wire edge. 

description of the power dependence. Particularly near 
the wire edge, higher pump powers were required exper- 
imentally to observe the same effects. Our energy re- 
laxation model assumes that the relaxation rate is pro- 
portional to the condensed polariton population and the 
energy difference in a relaxation process. In reality, we 
expect the rate to also depend on the total exciton popu- 
lation and not to be such a simple function of the energy 
difference. The study of these issues is beyond the scope 
of our work; it would require further experimental char- 
acterization and a microscopic model of the polariton- 
reservoir kinetics. Although we also neglected exciton 
diffusion, intrinsic reservoir dynamics and nonlinear loss 
mechanisms, we believe that our model captures the most 
important elements of the physics involved and repro- 



duces the experimental character. 



IV. CONCLUSIONS 

We have presented a theoretical formalism that models 
both the condensation process, as polaritons relax from 
high-energy hot-exciton states into the condensed polari- 
ton field, as well as further energy relaxation taking place 
within the condensate itself as polaritons interact with 
excitons and phonons. The theory goes beyond the stan- 
dard Gross-Pitaevskii theories, which do not describe the 
condensation process or energy relaxation, and the Boltz- 
mann rate equations, which do not allow a description of 
coherence, interference or the structure of modes in con- 
fined systems. 

We applied the formalism to the recent experiment of 
Wertai et al. This experiment uncovered a number of 
remarkable effects in the physics of nonequilibrium Bosc- 
Einstcin condensation in confined systems. The presence 
of optically induced potentials changes the modal distri- 
bution in the system. Condensation initiates strongest 
locally, at the position of the optical pumping. However, 
the dispersion of polaritons and their further relaxation 
allows lower energy modes to host the condensate in a 
spatially dependent way. Our model reproduces well the 
experimental findings. In particular it explains the puz- 
zling feature of discrete energy modes of the condensate, 
observed experimentally. The excitation of this mode 
ladder is naturally attributed to the energy relaxation 
mechanism within the condensate. A selection rule in 
the energy relaxation steps is determined by the balance 
between gain and loss for the final energy mode of the 
relaxation process. In the case of excitation at the wire 
edge, the narrow laser-induced trap creates a clear set of 
harmonic modes. 
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Quantum Photonics (NCCR QP), research instrument of 
the Swiss National Science Foundation (SNSF). 
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